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Abstract — In this paper, we study the performance of an 
amplify-and-forward (AF) based analog network coding (ANC) 
relay scheme in a multi-hop wireless network under individual 
power constraints. In the first part, a unicast scenario is con- 
sidered. The problem of finding the maximum achievable rate 
is formulated as an optimization problem. Rather than solving 
this non-concave maximization problem, we derive upper and 
lower bounds for the optimal rate. A cut-set like upper bound 
is obtained in a closed form for a layered relay network. A 
pseudo-optimal AF scheme is developed for a two-hop parallel 
network, which is different from the conventional scheme with 
all amplification gains chosen as the maximum possible values. 
The conditions under which either the novel scheme or the 
conventional one achieves a rate within half a bit of the upper 
bound are found. Then we provide an AF-based multi-hop ANC 
scheme with the two schemes for a layered relay network. It 
is demonstrated that the lower bound of the optimal rate can 
asymptotically achieve the upper bound when the network is in 
the generalized high-SNR regime. In the second part, the optimal 
rate region for a two-hop multiple access channel (MAC) via AF 
relays is investigated. In a similar manner, we first derive an 
outer bound for it and then focus on designing low complexity 
AF-based ANC schemes for different scenarios. Several examples 
are given and the numerical results indicate that the achievable 
rate region of the ANC schemes can perform close to the outer 
bound. 

Index Terms — Amplify-and-forward (AF), analog network cod- 
ing (ANC), wireless networks, multiple access. 



I. Introduction 

SINCE the introduction of the amplify-and-forward (AF) 
relay scheme, it has been studied in the context of coop- 
erative communication [l]-[4]. It is an interesting technique 
from the practical standpoint because the complexity and cost 
of relaying, always an issue in designing cooperative networks, 
is minimal for AF relay networks. As the simplest coding 
scheme, AF is also used to estimate the network capacity of 
relay networks. Obviously, the achievable rate of AF scheme 
can be viewed as a lower bound to the network capacity. In 
addition to its simplicity, AF is known to be the optimal relay 
strategy in many interesting cases [5]-[7]. Gastpar and Vetterli 
[7] have shown that for a two-hop network with AF relays, 
the cut-set bound of the network capacity can be achieved in 
the limit of a large number of relays. 
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Network coding [9] is a novel and promising design 
paradigm for wired communication networks. As opposite 
to the conventional routing operation, network coding allows 
the intermediate relays processing the received packets to 
reduce the amount of transmissions and thus improves the 
total throughput of the network. Li et al. [8] have shown that 
linear network coding can achieve the multicast capacity [9] 
in a noiseless network. This result indicates that to send out 
a linear combination of the incoming packets at each node is 
sufficient to obtain the optimal capacity performance. Linear 
network coding has also been studied in [10] from an algebraic 
perspective. Each destination node effectively obtains source 
information multiplied by a transfer matrix consisting of global 
encoding kernels on the incoming edges of it, and can recover 
the original data provided that the matrix is invertible. 

Applying the principle of network coding to wireless com- 
munication networks has recently received tremendous atten- 
tion from the research community. AF relay scheme allows 
one to exploit the broadcast nature of the wireless medium 
and to introduce the concept of network coding into physical 
layer. Katti et al. [11] have studied an AF-based analog 
network coding (ANC) scheme. As opposite to the traditional 
approach of wireless communications, analog network cod- 
ing fully use the interference rather than avoiding it. This 
technique significantly improves the network throughput in 
many scenarios. Since then, many works have been focused 
on the design of ANC relay schemes both for one-way and 
two-way relay channels [12]-[17]. In [15] Marie et al. have 
studied a multi-hop ANC scheme for a layered relay network 
under the individual power constraint. With each relay node 
amplifying the received signal to the maximum possible value, 
the achievable rate is shown to approach the network capacity 
in the high-SNR regime. Later such result is extended to the 
generalized high-SNR regime as defined in [16]. The results in 
[16] have shown the effectiveness of ANC in such generalized 
scenario. Recently, Agnihotri et al. [17] propose an iterative 
algorithm to obtain the optimal ANC scheme for a layered 
relay network with "Equal Channel Gains between Adjacent 
Layers (ECGAL)" property in general SNR regime. 

To employ ANC in a wireless network with AF relays, we 
are especially interested in deriving the optimal achievable 
rate. There have been many results for this problem [18]- 
[20]. For a two-hop parallel relay network, Marie and Yates 
[18] have found the optimal AF relay scheme in closed-form 
along with the maximum achievable rate under the sum power 
constraint. A more general result is obtained by Gomadam 
and Jafar [19]. A case when all the relay nodes introduce the 
correlated Gaussian noises is considered. They have found the 



optimal AF relay scheme for this scenario and investigated the 
influence of the correlation between noises on the end-to-end 
performance. However, the optimal AF scheme for a general 
network where both the topology and the operation regime 
of it have no limitation, is still unknown, even in the case of 
a two-hop parallel relay network under the individual power 
constraints. Agnihotri et al. [20] have provided a framework 
to compute the maximum achievable rate with AF schemes 
for a class of general wireless relay networks, which casts 
the problem as an optimization problem. The similar idea 
also appears in an independent work [16]. Unfortunately, the 
optimization problem in general case is hard to be solved. 

To employ AF-based ANC scheme in a wireless network, 
there are still many important and interesting problems re- 
maining to be solved. We list some of them as follows. 

• Problem 1: The cut-set bound is usually used to justify 
the performance of different relay schemes. However, it 
is not a tight upper bound for the AF achievable rate in 
most cases because the relay nodes are only allowed to 
do linear operations. To derive a tight upper bound for the 
achievable rate of a network with AF relays is an issue 
of great importance since the optimal rate is always hard 
to be obtained. 

• Problem 2: Is it always optimal for relay nodes to 
amplify the received signals to the maximum possible 
values under individual power constraint? Intuitively, for 
different scenarios there may exist different AF schemes 
that have a better end-to-end throughput performance than 
any others. Then how to characterize them? 

• Problem 3: A multi-user scenario is also an interesting 
topic [21]. For a two-hop MAC with AF relays, how 
many AF schemes shall we use to obtain the entire 
achievable rate region? Obviously, it is infeasible from 
the practical standpoint for the relay nodes to store all the 
AF schemes. Hence, it is worth investigating the tradeoff 
between the performance and the complexity. 

The main results of this paper give partial answers to the 
above questions and are summarized as follows: 

• In this paper, we pursue two related objectives. The first 
one is to investigate the optimal AF-based ANC scheme 
for a multi-hop unicast relay network. Assuming each 
relay node has a transmitting power constraint, we derive 
an upper bound of the optimal achievable ANC rate in a 
closed form for a layered relay network. The idea behind 
our method is similar to the technique used to derive the 
cut-set bound. Thus we call it a cut-set like upper bound. 
A novel AF scheme for a two-hop parallel relay network 
is proposed, which is different from the conventional 
scheme with all relay nodes amplifying the received 
signal to the maximum possible value. We determine 
the different conditions under which either the novel 
scheme or the conventional one has a better performance 
to approach within 1/2 bit to the upper bound and thus 
the optimal rate. Based on such observations, a mixed 
AF-based multi-hop ANC scheme is proposed. When the 
network is in the generalized high-SNR regime [16], we 
demonstrate that the lower bound of the achievable rate 



can asymptotically achieve the upper bound. Furthermore, 

we point out that the result obtained in [15] can be viewed 

as a special case of the result obtained in our paper 

• The second objective of this work is to extend the ANC 

scheme to the multiuser case. S. A. Jafar et al. [21] have 

shown the optimal rate region of the two-hop parallel 

AF network under a sum power constraint. We find 

that to assume the relay nodes under individual power 

constraints makes the problem more complicated. It is 

observed that for a specific AF scheme, the network 

reduces to a conventional MAC. Each Gaussian capacity 

in the rate set of the MAC is equivalent to the capacity 

of a two-hop parallel unicast network. With the similar 

idea, we first derive two outer bounds for the optimal rate 

region. The structures of them are fully characterized. 

Then a dynamic ANC scheme is proposed. According to 

different network settings, the number of the AF schemes 

stored in the relay nodes varies and with finitely many 

schemes a time-sharing inner bound is obtained. The 

tradeoff between the storage space and the performance 

is decided for the practical purpose. Finally, we illustrate 

three ANC schemes. The numerical result shows that 

the inner bound of each ANC scheme is close to the 

respective outer bound in different scenarios. 

Notation: Scalars are denoted by lower-case letters, e.g., 

X, and bold-face lower-case letters are used for vectors, e.g., 

X, and bold-face upper-case letters for matrices, e.g., X. 

In addition, tracepC), |X|, X-^ and X^^ denote the trace, 

determinant, transpose and inverse matrix of X, respectively, 

and diag{xi, • • • , x„) denotes a block-diagonal square matrix 

with xi, • • • , x„ as the diagonal elements. X(i, j) denotes the 

{i,j)-th element of X, and rankpC) denotes the rank of X. 

||x|| denotes the Euclidean norm of a vector x. E\-] is the 

expectation operation. log{-) denotes the logarithm in the base 

2 and ln{-) denotes the natural logarithm. We use conv (•) and 

col {■) to denote the convex hull and the closure of a set. 

II. Network Model 

In the first part of this paper, we analyze a multi-hop 
Gaussian relay network with a single source-destination pair 
{S,D), which is represented by a directed graph Q — {V,£) 
depicted in Fig. 1. Assume each non-source node k introduces 
a Gaussian noise, all the relay nodes work in a full-duplex 
mode, and there is no circle path in the network. At instant 
n, the channel output at node fc is a linear combination of 
all noisy signals transmitted from its upstream nodes and the 
Gaussian noise introduced by itself and can be expressed as 

Vk [n] = ^ hj^kXj [n] + Zk [n], (1) 

jev(fc) 

where hj,k denotes the channel gain from node j to node 
k, V{k) represents neighboring nodes of node k which have 
one-hop path to node fc, Xj [n] is the transmitted signal 
at node j, and Zk [n] is a sequence of independently and 
identically distributed (i.i.d.) Gaussian noises with zero mean 
and variance 1. All the channel gains are supposed to be fixed 
positive real-valued constants known through the network for 
the scope of the present paper. 




Fig. 1. Wireless relay network with single source-destination pair 



A natural assumption is that there exists a transmitting 
power constraint Pf, ^ at each node k such that 



E [4] < Pi!'. 



(2) 



Each network node performs analog network coding via 
amplify-and-forward relay scheme. Assume the relay nodes 
operate instantaneously as in [12], that is, the relay nodes 
amplify and forward their input signals without delay. 



Xk [n] = Pkyk [n] 



(3) 



where the amplification gain /3j. is chosen such that the power 
constraint (2) is satisfied. Actually the possible system insta- 
bility resulting from this assumption is avoided by a "buffering 
and subtracting" technique as observed in [20]. So, with this 
assumption, the time index n then can be omitted in the sequel 
for the sake of brevity. Through a relaying path from node j to 
node k, the signal Xj is multiplied by the amplification gains 
and the channel gains along this directed path. Since the multi- 
hop ANC scheme takes more of the network coding approach, 
the corresponding coding coefficients are defined as follows. 
Definition 1 (Local Encoding Coefficient): Let e{i,j) and 
e{j, k) be the input and output channels of node j. We call 



e{j,k) — Pj'''e{j,k) 



(4) 



the local encoding coefficient of pair {e{i,j),e{j,k)),i G 
V{j), where, in particular, /3s = 1- 

Definition 2 (Global Encoding Coefficient): We call 



f3,k 



n 






(5) 



the global encoding coefficient from node j to node k, where 
£{j, k) represents the set of channels appearing in a relay path 
between the two nodes. 

The global encoding coefficient defined here can be inter- 
preted as the equivalent path gain between node j and k. The 
relationship between the local and global coefficients is given 
as follows. 

fs\k = 2^ 0^e(j,k)fs.J (6) 

iev(fe) 

III. Upper Bound to ANC Rate 

A. Sufficient Condition of Power Constraint 

We first derive the expression of the achievable rate for 
a unicast network with a specific ANC scheme. Then the 
problem of finding the maximum achievable rate is formulated 



as an optimization problem. In this subsection, we focus 
on capturing the feasible domain of each amplification gain 
according to the power constraint at each relay node. 

According to the assumption that each non-source node 
introduces a Gaussian noise, each of them can be viewed 
as a source with the noise as the input. Therefore, by the 
linear operation of AF relay scheme, the original network 
can be considered as a linear combination of several unicast 
subnetworks and the source of each subnetwork is either the 
source node S or the non-source node. Therefore, with the 
global encoding coefficient (5), the received signal at the 
destination is expressed as 



yo = fs,DXs 



•^S,eq 



{j} 

^eq-) 



(7) 



where xs.eq — fs,DXs represents the equivalent signal trans- 
mitted from the source node, and z^q — X) fjDZj + zd 

{3}' 
denotes the equivalent noise received at the destination node 

where the summation is over all relay nodes {j}. By the 
assumption that all the Gaussian noises introduced at the non- 
source nodes are independent, the sum of them is also drawn 
according to Gaussian distribution. Hence, from the RHS of 
the second equality above, the relay network employing the 
ANC relaying scheme can be considered as a point-to-point 
Gaussian channel. As the well-known result, the source node 
adopts the Gaussian codebook. The random variable (rv.) 
Xs used to generate the codewords is drawn according to 
N [0, Ps] , where N [O, cr^] denotes normal distribution with 
zero mean and variance cr^ . A sequence of codes containing 
|"2n-R"| codewords of length n is proposed and it is shown that 
the error probability goes to zero as n ^ 00. 

From (4), (5) and (7), both the equivalent signal and noise 
are related to the amplification gains chosen in the ANC 
scheme. For the equivalent Gaussian channel, we focus on 
the choice of the amplification gains to maximize the SNR 
value at the destination node, and thus the transmission rate 
R. Certainly, this problem can be formulated as a standard 
optimization problem with the object function being the SNR 
received at the destination node which is shown as follows. 



SNR{{Pk,k(^V\{S,D)}) = 



HdPs 



L^ Jj,D 

3eV\{S,D) 



fln + r 



(8) 



which is subject to the power constraints at the relay nodes. 
Hence, to solve this optimization problem, it is necessary to 
convert the power constrains into the constraints with respect 
to the amplification gains first. However, for a multi-hop relay 
network, the selection of the amplification gain at one node 
also relies on the selections at its upstream nodes. Therefore, 
we find it hard to give an equivalent independent constraint 
for each amplification gain. To simplify the problem, we first 
provide such an independent constraint. Then we show it 
suffice for the power constraint. Fortunately, it will be shown 
that under such constraint, a good multi-hop ANC scheme can 
be proposed in a special SNR regime defined in the sequel. 
Several parameters are defined first as observed in [15]. 



I l+l 



L-\ 




Fig. 2. Layered relay network with single source-destination pair 

Definition 3: When each node j transmits with P- ^ given 
in (2), the power received at node k is upper bounded by 

2 



Pi 




R,k 



and the reciprocal of Pr ^ is represented by 

1 



P; 



R.,k 



(9) 



(10) 



Note that if the network is noiseless, then the received signal 
power at node k approaches (9) when all the neighboring 
nodes of it transmit the received signals at the upper bounds 
of the power constraints (2). Then we extend the sufficient 
condition of the power constraint in [15] to a general network 
model described in section II by the following lemma. 

Lemma 1: In a multi-hop Gaussian relay network, it is 
sufficient for each relay node k to choose the amplification 
gain 

pUp 



fit< 



(l + 4)Pfl,, 



(^. 



(11) 



such that the power constraint (2) is satisfied. 

Proof: The proof is given in the appendix A. ■ 

Note that by Definition 3, the value PR^k and its reciprocal 

5k are both constants determined by the network settings. 

Therefore, the constraints for the amplification gains proposed 

in the above lemma are independent of each other. 

B. Upper Bound to ANC Rate 

In this subsection, we first show a layered relay network 
depicted in Fig. 2. The source node is assumed to be at layer 
and the destination node at layer L. The number of relay 
nodes at layer / is denoted by rij. The transmitted signal vector 
at layer I is denoted by x; and the received signal vector at 
layer / + 1 can then be obtained as 



yi+l = H;X; +Z;+i, 



(12) 



where yi+i and Z(+i are the received signal vector and the 
noise vector at layer / + 1 respectively, and H; denotes the 
channel matrix between these two layers with the element 
H; (j, k) representing the channel gain from node j at layer I 
to node k at layer l + l. Especially, we use ho to denote the 
broadcast channel between the source node and the nodes in 
the first layer and hi_i to denote the multiple access channel 
between the nodes in the L — 1th layer to the destination node. 



Specifically, the global encoding coefficients in such net- 
work are given as follows. 



fs,D — i^L-i^L-1 ■ ■ ■ HiBiho, 



(13) 



/: 



J,D 



h^iBL-i • • • hi,,/3„ j e A, / = 1, 2 • • • i - 1, (14) 



where h; j represents the jth column-vector in H;, and B; = 

diag {/3i ■ ■ ■ /3j ■ ■ ■ /?„, } ,j G Ci where £; denotes the node 
set of layer I. From (7), (13) and (14), the signal received at 
the destination node is given by 

VD = h|^_iBL_i • • • HiBihoXs 

L-l 

+ ^hLiBL-i---H,B,z,+zz5, (15) 
1=1 

and the SNR function at the destination node (8) results to be 

SNR{{lik,kGCuC2---CL-i}) 



E 



(hLiB 



L-l • • • 



HiBihoXs 



E 



E^hLiBi-i-'-HiBiZz 
1=1 



(16) 



E[zl] 



To have an upper bound to the ANC rate, we recall the 
technique of obtaining the cut-set bound of capacity of a noisy 
network. One may first assume all channels in the network are 
noiseless except at one cut and then have an upper bound. By 
taking the minimum upper bounds over all cuts, the cut-set 
bound is derived. By applying the same idea to ANC in layered 
networks, an analysis of an ideal layered relay network is made 
first. We first fix a layer la and assume the Gaussian noises 
are only introduced at nodes on this layer and other parts in 
the network are noiseless. By optimizing the SNR function 
received at the destination node of such network which is 
formulated by 



SNRi„{{l3k.k(^Ci,C2---CL-i}) 



E 



{\i^_-^Bl-i ■ ■ ■'BiYiQXs 



E (h|;_iBi_i...H,„B,„z,J 



(17) 



an upper bound of the optimal ANC rate Ropt is derived. Then 
by traversing ?o = 1,2 •••L, a collection of upper bounds 
is obtained. To choose the minimum of them, a better upper 
bound is derived. We draw the conclusion in the following 
theorem. 

Theorem 1 (Upper Bound to ANC Rate): An upper bound 
of the ANC rate of a Gaussian layered relay network with L 
layers under the individual transmitting power constraints is 
given as follows. 

R"""^, min ilog(l + PL,Pi^,;o), 



where T'n^i 



J e Ao- 



VPr1---VPrJ---Vp^o_ 

Proof: The proof is given in the appendix B. 

We refer to the optimal AF scheme at layer Iq for the 

ideal network with only the nodes in layer Iq introducing 

the noises as a "pseudo-optimal amplification scheme". Note 

that the idea behind our upper bound is similar to the idea of 



cut-set bound, we call this upper bound a cut-set like upper 
bound. As we know, the cut-set bound is an upper bound to 
the network capacity, where the relay nodes are allowed to do 
any operations, including decoding. So we intuitively feel that 
the cut-set bound may not be tight for the optimal achievable 
rate of a network with AF relays in general case. To generalize 
the result of Theorem 1, we hope to find out a more general 
upper bound result for a non-layered relay network with AF 
relays. By taking a cut set of the network, we use the notation 
J to denote the nodes on the output sides of the channels 
in such cut set. If assume only the nodes j's in J introduce 
the Gaussian noises, the corresponding SNR function can be 
expressed as. 



SNRj{{^^,,k^V\{S,D)}) 



flpPs 



(18) 



leJ 



0,D 



Using the same argument as we derive Theorem 1, by opti- 
mizing SNRj, an upper bound of Ropt is derived. Then the 
upper bound of ANC rate in general network is given as 



R^P — min max - log ( 1 

J l3k,k£V\(S,D) 2 



SNR 



J) 



(19) 



However, we are unable to give a closed-form of the optimal 
solution for SNRj in general case. 

Remark 1: The upper bound of the ANC rate in [15] can be 
regarded as a corollary of Theorem 1. Let us revisit the result 
obtained in [15] in an alternative approach proposed in our 
paper. Since in the scenario considered in paper [15], all the 
relay nodes have sufficiently large transmitting powers except 
the ones at layer L — 1, the upper bound in Theorem 1 shows 
that the minimum is taken when /q = L. The SNR function of 
an ideal network with Iq — L , i.e., only the destination node 
introduces the Gaussian noise turns to be 



SNRL{l3k,keCi,C2---CL-i} 



E 



(hLiB 



L-l • • 



•HiBihoXs 



E[zl] 



(20) 



then the corresponding upper bound of the achievable ANC 
rate is 

]-\og{l + Pr,d)- (21) 



R^P 



Furthermore, from (20) we can see that the noise power 
received at the destination is independent of amplification 
gains in this case. Hence, the larger the transmitting powers 
are, the better the performance of the ANC scheme. With 
the power constraint, a practical ANC scheme is proposed as 
observed in [15]. 



oUp 



/3fc = /3fc^fceA,/ = l,2,--- ,i-l. 



(22) 



Remark 2: In [7], an upper bound to capacity of the relay 
network is computed using a weaker corollary of cut-set 
theorem [25, Theorem 14.10.1] and the capacity formula for 
Gaussian vector channels with fixed transfer function [24]. 
Applying the similar idea used to verify the corollary, we can 
also compute the upper bound to capacity of the layered relay 



network, which can be formulated as an optimization problem 
as follows. 



mm 



max /(x,„_i;y,„) 



subject to 



E[x)]<PfP,3^Ci,^, 



Let us relax this power constraint to 



E ^H]^ E p. 



Up 



(23) 



(24) 



(25) 



ie^io-i 



je^io-i 



Note that for given ^o, (23) is formulated due to the assumption 
that the nodes in layer Zq — 1 are the multiple transmitting 
antennas of the source node and the nodes in layer Iq are the 
multiple receiving antennas of the destination node. Hence, 
the transmitting antennas may have any cooperation to encode 
the message and the multiple receiving antennas may have 
any cooperation to decode it. Consequently, the corresponding 
capacity is an upper bound to the network capacity. The 
solution of the problem defined by (23) and (24) cannot be 
larger than the solution of the problem defined by (23) and 
(25). For a specific l^, the latter can be evaluated by using the 
result in [24]. 

Ci„ = ^max -log (det (/„,^ + H,„_iKx,„_iH?JJ_i)) , 

(26) 
where Kxj _^ ~ E [x/j,_ixf _j^] and the maximum is subject 
to trace (Kxij,_i) < E Pf^- Note that the optimal 

solution is obtained by water- filling. To investigate how far 
between the cut-set bound and the upper bound derived in 
Theorem 1, several special cases are taken into account. 

When the minimum in Theorem 1 is obtained by taking 
Iq = 1, i.e., the first hop is on the bottleneck in the data 
transfer, we have the broadcast cut-set bound as the upper 
bound for capacity. 



Cbc = \ log (1 + hJhoPs) 



(27) 



At this time, the upper bound obtained in Theorem 1 is 

1 



i?^P--log(l + h^hoPs) 



(28) 



which is equal to the former. Hence, in case when i?'^^ is 
achieved, then the network capacity can be achieved. 

When the minimum in Theorem 1 is obtained by taking 
Iq — L, i.e., the last hop is on the bottleneck of the data 
transfer, we have the multiple access cut-set bound as the upper 
bound for capacity. 



Cmac 



1 



■log i + h£_ihi_i 



2log(l 



iL-ihi- 





/jjj = 1 ^^^^\ h^^ = 1 




1^20= a 



D D 

Fig. 3. Two-hop parallel relay network with single source-destination pair. Fig. 4. 3-layer network with 2 nodes at layer 2. 



where P^p == 


V^i"^---v/^"'--V<- 


. At this time, 


the upper bound obtained in Theorem 1 is 




R^P = ^log{l + PnM) 






= liog(i + (hLiP^^)^) 


(30) 



By Schwarz inequality, R^p < Cmac- The equality holds 
when V^P — ch^-i for some constant c. 

IV. Lower Bound to ANC Rate 

We expect to obtain a good AF-based ANC scheme by 
solving an optimization problem in this section. The SNR 
function given in (16) is the object function, which is subject 
to the constraint given in Lemma 1. Note that Lemma 1 
only provides a sufficient condition for the power constraint. 
Therefore, the optimal value of this optimization problem can 
only be served as a lower bound to the optimal ANC rate. 



max SNR{{(ik,k^Ci,C2---CL-i}) 



s.t. 



Pl< 



Pk 



PR.k + 1 



,k e Ci,C2- ■ ■ Cl-1 



(31) 



As shown in [18], a similar optimization problem has a 
simple solution under the sum power constraint for a two-hop 
relay network as depicted in Fig. 3. However, the optimization 
problem given above is essential a fractional programming 
[22] which is a hard problem. Even for a two-hop network, 
the above problem is a fractional programming [23] where 
both the enumerator and denominator are convex functions. 
There is no such easy approach to solve the problem so far as 
we know. 

Before investigating the lower bound to the optimal ANC 
rate for the layered relay network with arbitrary layers, we 
first examine three different AF schemes for two single-source 
single-destination two-hop parallel relay networks as shown in 
Fig. 4, and Fig. 6. Then we capture the characterizations of 
these schemes. 

• Scheme 1: As observed in most conventional amplify- 
and-forward relaying schemes, all the relay nodes trans- 
mit the received signal at the upper bounds of the power 
constraints. Hence, the amplification gains are 



Pk 



aV.^ 



l,2---n. 



Scheme 2: Motivated by the scheme used in layer Iq when 
we obtain Theorem 1, we employ the pseudo-optimal 
amplification scheme for this two-hop network, that is 

h = c- — ^JPR.k.k == 1,2-- -n, 
rikD 



where c = min < /3^ 



Vv 



i-kD 



l,2---n 



Scheme 3: The destination node chooses one relay node 
to transmit the received noisy signal which is referred to 
as AF with selection. To be specific, only the relay node, 
with the assistant of which the destination receives the 
highest SNR value, amplifies and forwards its received 
noisy signal with the upper bound of the power constraint. 
Find 

fco = arg ma.x( hskhkoPk^ 



then choose 



fe=l,2---n 



Mp 



h, = Pko and /3fc - 0, fc ^ fco. 
For the network parameters given in Fig. 4, we compute 
the corresponding amplification gains and the SNR values 
received at the destination node for the above three schemes. 
Assume each Gaussian noise introduced by the non-source 
node is i.i.d. with unit variance. Later, we will explain that 
the assumption does not lose the generality. 

• Scheme 1: 



Scheme 2: 
For < a < 1, 
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V2 



and SNR = 
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and for 1 < a. 
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Scheme 3: 
For < a < 1, 
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and for 1 < a. 
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Fig. 5. Achievable rates of different AF scliemes. 



Then we compute the upper bound of the achievable rate 
for this network. By Theorem 1, 



R^P = 



{rX\ 



R^,^}, 



where i?^^ = 0.5 log 3 and R^^ = 0.5 log [l + (a' 

The result shows that if < a < 1.3028, then R^ ^ is greater 
than i?2 ^, which means the bottleneck of data transfer is on 
the second hop, and if 1.3028 < a, then R^ ^ is no larger than 
i?2 ^, which means the bottleneck of data transfer is on the 
first hop. 

In Fig. 5, the numerical result of the achievable rate for 
each scheme is given. From Fig. 5, we observe that when 
< a < 2, the achievable rates of scheme 1 and 2 are almost 
the same, while when 2 < a, the achievable rate of scheme 
2 outperforms that of scheme 1. We are interested in how to 
select the AF scheme in different cases. Since in low SNR 
scenario, it is not expected to perform AF as the relaying 
scheme, to investigate the criterion, we consider the case when 
a is sufficiently large, i.e., a > 2. The total noise received at 
the destination node has two portions, one is the propagation 
noise from relay nodes denoted by wd, which is referred to 
as the upper layer noise hereafter, and the other is the noise 
introduced by the destination node itself which is referred to 
as the lower layer noise later We compute the propagation 
noise power both for scheme 1 and 2. In scheme 1, 



and in scheme 2, 



E[i 



+ a% 



'D\ 



Since the lower layer noise has unit power, we observe that 
in both cases, the propagation noises dominate the total noise 
received at the destination node. Therefore, the AF scheme 
proposed in our paper is expected to be more effective in 
dealing with the upper layer noises. Note that the dominance 
only depends on the ratio of the variances of the upper layer 
and lower layer noises, hence assuming all the Gaussian noises 
have unit variance is without loss of generality. 
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Fig. 6. 3-layer network with 2 nodes at layer 2. 
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Fig. 7. Achievable rates of different AF schemes. 



For the network parameters given in Fig. 6, we compute 
the corresponding amplification gains with the above three 
schemes. 

• Scheme 1: 



VT+o^ vTT^ 



and 



SNR = 



1 



• Scheme 2: 
For < a < 1, 

Pi = ,^ Q . /32 = , and 



a^/l + a? 



Vm? " ^/l + ^ 
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2(g2 + 1) 



and for 1 < g, 
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• Scheme 3: 
For < g < 1, 
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and for 1 < g. 
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Then we compute the upper bound of the achievable rate 
for this network. By Theorem 1, 



i?^P=min{i?fP,i?2^P}, 



where i?f ^ == 0.5 log (l + a^ + a^) and R^^^ = 0.5 log 5. The 
result shows that if < a < 1.2496, then R^^ is less than 
i?2 ^, which means the bottleneck of data transfer is on the 
first hop, and if 1.2496 < a, then R^ ^ is no less than R2 ^, 
which means the bottleneck of data transfer is on the second 
hop. 

In Fig. 7, the numerical result of the achievable rate for 
each scheme is given. From Fig. 7, we observe that when 

< a < 1, the achievable rates of scheme 1 and 2 are almost 
the same while when 1 < a, the achievable rate of scheme 

1 outperforms that of scheme 2. We are still interested in 
how to select the AF scheme in different cases. As explained 
previously, we consider the case when a is sufficiently large, 
i.e., a > 1. We compute the propagation noise power both for 
scheme 1 and 2. In scheme 1, 



E w'n = 



1 



1 



and in scheme 2, 



1 



^ [--] - iT-a 



i4 a2(l + a4)' 

We observe that in both cases, the lower layer noise dominates 
the total noise. Therefore, from the numerical result, the 
scheme 1 is expected to be more effective in dealing with 
the lower layer noise. 

Now let us summarize the results drawn from the above 
analysis. The first observation is that the more relay nodes 
assist the transmission the better performance we can obtain. 
The reason is that the source signals retransmitted by the 
relay nodes are all related; however, the noise signals are 
independent of each other Therefore, from the aspect of 
cooperation communications, the more copies of the original 
signals received at the destination node the lager diversity 
gain will be obtained. Second, if the relay nodes always 
amplify and forward the received signal to the maximum 
possible value, then it may result in suboptimal end-to-end 
throughput. This observation leads to an interesting problem 
that under what conditions the two schemes will have a better 
performance? From the above two examples, we conjecture 
that scheme 1 is more effective in managing lower layer noises 
and inversely, scheme 2 is more valid in dealing with upper 
layer noises. Finally, the first two AF schemes mentioned 
above can constitute a mixed ANC scheme which can be 
applied to a layered relay network with arbitrary layers. The 
relay nodes at each layer may select either of the conventional 
scheme (scheme 1) or the new scheme (scheme 2) proposed in 
our paper. However, in a layered relay network, the selection of 
the scheme for one layer also depends on the schemes applied 
by the other layers. Therefore, it is complicated to provide a 
localized criterion to select which scheme at each layer 

Let us consider the following AF-based ANC scheme mo- 
tivated by the optimal solution in the ideal network analyzed 
in Theorem 1 . With such ANC scheme, a lower bound to the 



optimal ANC rate is derived. To simplify the computation of 
the lower bound, we first define 



5{) = max 5j. 

jeCi,--- ,Cig-i,Ci„^i,--- ,Cl 



(32) 



The amplification gain at node k,k E £i,£2 ■ ■ ■ ■Cig-i is 
chosen as 



/3fc- 



P, 



Up 



< 



p; 



Up 






and at node k,k E Cio+i • • • 'Cl-i is chosen as some positive 
constant such that 



/3fc — const. < f3f, ^ . 



(34) 



Note that the amplification vector selected at layer Iq depends 
on these factors and moreover, no matter what feasible gains 
are selected will not affect the resulting upper bound. The 
amplification vector at layer Iq is chosen as 



Ao — ciG Pfl,/o, 



(35) 



where /3;„ and G are defined in the proof of Theorem 1, and 
ci can be determined as 



Cl 



min < /?j, ^ 



9k 



VP^ 



,ke Ci„ 



(36) 



Note that, in the above ANC scheme, both the conventional 
and the pseudo-optimal amplification scheme are adopted. 
Although it is not expected to obtain the optimal ANC rate 
through a fixed scheme, the corresponding achievable rate of 
such ANC scheme can be served as a lower bound. 

Theorem 2 (Lower Bound to ANC Rate): A lower bound 
of the ANC rate of a Gaussian layered relay network with L 
layers under the transmitting power constraints is as follows. 
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Proof: The proof is given in the appendix C. ■ 

Note that C2 is a bounded positive constant depending on 
the selection of amplification gains in nodes at layer ?o + 1 
to i — 1. Finally, we should emphasize that the lower bound 
obtained in the above theorem is a universal one since it is 
not required that (5o should be a small positive value. 

V. Applications of AF-based ANC scheme 

To the best of our knowledge, AF relay design for a layered 
relay network has been the focus of mush recent research. 
Thus by taking pseudo-optimal amplification scheme at certain 
layers, we have an alternative multi-hop ANC scheme as 
opposite to the scheme with all relay nodes amplifying the 



received signals to the maximum [15]. A problem is whether 
pseudo-optimal amplification scheme can make things better 
than conventional AF scheme. We turn to two particular 
classes of layered relay networks and compare the resulting 
lower bound of AF rates to our upper bound. 

A. Application to Two-Hop Network 

We extend the results derived from the previous examples 
to a two-hop parallel relay network with n relay nodes. Note 
that scheme 2 can be viewed as a pseudo-optimal amplification 
scheme for such network. Denote by B*^^' and B*^^^ the 



„(i) 



(2) 



two schemes and by w}^ and w)j the upper layer noises 
corresponding to these two schemes respectively. Hence, we 
have 



E 



(w 



Wx2 
D J 



hfB(*)BWhi,i = l,2, 



(37) 



and the lower layer noise introduced by the destination node 
itself has, as assumed, a unit power 



E zi, =1. 



It is easy to verify that E 



(w 



(1)^2 



(38) 
is no less than 



(-g^)^ 



E 

Definition 4: We call the upper layer noise dominates the 
total noises received at the destination in a two-hop relay 
network with AF relays when 



E 



i<'r 



> E 



(2) 



«') 



> 1. 



Reversely, we call the lower layer noise dominates the total 
noises when 
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(2)^2 

D 



i<^) 



< E 



(w 



(1)n2 
D ) 



< 1. 



As we mentioned before, the individual power constraints 
make the problem of deriving the exact optimal AF rate 
complicated to be handled. However, we still obtain some 
interesting results on the performance of achievable AF rate 
shown as follows. 

Theorem 3: For a two-hop relay network with AF relays 
under individual power constraints, there exist AF schemes 
accodrding to different senarios when either the upper layer 
noise or the lower layer noise dominates the total noises 
recived at the destination node such that the gaps between 
the the corresponding achievable rates i?'s and the optimal 
AF rate Ropt are at most 1/2 bit. 

Proof: For a two-hop relay network, we see that the 
achievable rate can be easily derived without employing the 
lower bound obtained in Theorem 2. For the first scenario, 
we compute the achievable rate with the pseudo-optimal 
amplification scheme, i.e., B*^^^. 



R 





hfB(2)hohjB(2)hi/s liw^a^f 



Ps 



hfB(2)hohjB(2)hiPs 



E 



I (2)\2 



where (a) follows from the assumption that E 
and (b) follows from Theorem 1 that 



(-g¥' 



(39) 
> 1 



2 \ E («.g')» I 



For the second scenario, we compute the achievable rate 
with the conventional scheme, i.e., B^^^ 




where (a) follows from the assumption that E {w\ 
and (b) follows from Theorem 1 that 



Ropt < 



^"2" = \ log (1 + hf B(i)hoh;5'B(i)hiP5 



Then we complete the proof. ■ 

Note that Theorem 3 rigorously interpret the intuition 
obtained from the examples in section IV. We are also in- 
terested in the case beyond these two scenarios, i.e., when 



E 



{^'^f 



> 1> E 



i<^? 



Are the performances of the 



two schemes in such general scenario also good enough? We 
draw the conclusion in the following theorem. 

Theorem 4: For a two-hop relay network with AF relays 
under individual power constraints, with either the pseudo- 
optimal amplification scheme or the conventional one, the 
gaps between the corresponding achievable rates _R's and 
the optimal AF rate Ropt are upper bounded by a constant 
when neither the upper layer noise nor the lower layer noise 
dominates the total received noises at the destination. 

Proof: To verify the statement, we compare achievable 
rates corresponding to scheme 1 and 2 with two upper bounds 



jt/p 



jC/p 



i?2 and R^ ^ respectively. The achievable rate with the 
conventional scheme B^^^ is as follows. 
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(42) 



where the last inequality results from the generalized Rayleigh 
quotient [27], and (c) follows from the same argument as (b) 
in (40). Similarly, the achievable rate with the pseudo-optimal 
amplification scheme B^^) j^ ^s follows. 
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where (a) follows from the assumption E 

follows from (42) and (c) follows from the same argument as 

(b) in (39). 

Therefore, we complete the proof. ■ 

B. Application to Generalized High-SNR Regime 

In this subsection, a special SNR regime is proposed, under 
which the ANC scheme proposed in section IV is shown to 
have a good performance. The scenario considered in this 
paper can be viewed as an extension of that discussed in [15]. 
In the similar manner as in [15], we define the generalized 
high-SNR regime as follows. 

Definition 5 (Generalized High-SNR Regime): We call the 
layered relay network is in the generalized high-SNR regime 
if 

So < S, (44) 

for some small positive S, where 6q is defined in (36). 

From the definition, it is assumed that the received SNR 
at each node k, k ^ Ci^ is large which is referred to as the 
high-SNR condition in the sequel. Therefore, the above SNR 
regime is an extension of the high-SNR regime defined in [15], 



where only the destination node does not satisfy the high- 
SNR condition. If we set lo equal to L, the generalized high- 
SNR regime is degraded to the high-SNR regime. From the 
perspective of the optimization problem, we will show that the 
AF-based ANC scheme proposed previously can be utilized 
as a good suboptimal solution for (31). It is expected that 
the achievable rate of such ANC scheme may asymptotically 
approach the upper bound of the ANC rate within a small gap 
and thus close to the optimal ANC rate. 

Two cases when the generalized high-SNR regime can be 
satisfied are considered. If assume except the nodes on layer 
Iq — 1, all the other nodes have sufficiently large transmitting 
powers and the transmitting powers of these nodes tend to 
infinity at the same rate as (5 ^ 0, the network is in the 
generalized high-SNR regime. Under these assumptions, we 
have PR.fe = const., k G £;„, C2 = const., and ci — > oo as 
(5^0. It is not complicated to verify that 
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The result implies that the achievable ANC rate asymp- 
tomatically approaches the upper bound. This analytical result 
can be justified intuitively as follows. The total noise received 
at the destination node is separated into three parts. We 
first make a comment on the part consisting of the noises 
introduced by the nodes on layer 1 to layer ^o ^ 1- As the 
source transmits the signal with large power, (if Iq ^ 1) and 
all the introduced Gaussian noises have finite powers, this part 
of noises seems to be small compared with the signal received 
at nodes e £;„ and also at the destination node. Then we 
discuss the second part consisting of the noises introduced by 
the nodes at layer Iq. As the generalized high-SNR regime 
indicated, the nodes on layer Iq — X have limited transmitting 
powers, each node on layer Iq receives a noisy signal with 
finite powers. Furthermore, as mentioned before, the propaga- 
tion noise is only a small fraction of it. Meanwhile, the noises 
introduced at layer ^o have powers on the same level as the 
received signal and hence this part of noises cannot be ignored 
at the destination node. Finally, we consider the last part of 
noises introduced by the nodes at layer /q + 1 to layer L. Since 
nodes at layer ^o possess large transmitting powers, the case is 
similar as the first part. Hence the introduced noises of these 
nodes are small compared with the transmitted signal from 
layer Iq, which implies that the propagated noises from layer Iq 
dominate the total noise received at the destination. Therefore, 
the corresponding achievable ANC rate may asymptotically 
approach the upper bound obtained from the ideal network 
under the assumption that only the nodes on layer Iq introduce 
the noises. 
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Fig. 8. 4-layer network with 2 nodes at layer 2. 



Next we discuss another case. Assume all the nodes both on 
layer Zq — 1 and ^o have finite transmitting powers, while any 
other nodes have sufficiently large transmitting powers. But we 
assume the number of nodes at the loth layer is sufficiently 
large which is similar to the assumption in [7]. Hence, there 
are indeed two modes the nodes satisfying the high-SNR 
condition. First, the nodes on layer 1,1 ~ 1 ■ ■ ■ Iq — I, lo+2 ■ ■ ■ L 
satisfy the high-SNR condition as the previous case. Second, 
the received signal powers at the nodes on layer ^o + l increases 
with the number of the nodes on layer Iq- Consequently, the 
generalized high-SNR regime can be satisfied. 

To emphasize on the difference between these two modes, 
let 

for some small S' > 0. Then we have. 
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The result also implies that the achievable ANC rate 
asymptomatically approaches the upper bound. Furthermore, 
we observe that with the number of the nodes in layer Zq 
increasing, the upper bound of ANC rate also increases and 
thus the optimal ANC rate. 

Finally, two examples are utilized to illustrate the perfor- 
mance of the ANC scheme in a 4-layer network for the 
generalized high-SNR regime with ^o = 2. The numeric results 
show that in both cases considered, the achievable ANC rate 
can asymptotically achieve the upper bound of the ANC rate. 
Example 1: 4-Layer Network with 2 relays at layer 2 
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Fig. 9. Acliievable ANC rate and upper bound of ANC capacity in 4-layer 
network. 
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Fig. 10. 4-layer network with n nodes at layer 2. 



A 4-layer relay network with 2 relay nodes on layer 2 
is depicted in Fig. 8. To satisfy the generalized high-SNR 
regime, assume relay nodes 3 an 4 possess sufficiently large 
transmitting powers. The lower bound of the achievable ANC 
rate is expected to increase with the transmitting powers of 
the source node. 

From Fig. 9 we observe the following: 

1) For Pr 3 and Pr4 constants, the achievable ANC rate 
asymptotically approaches the upper bound of ANC rate. The 
achievable ANC rate approaches the capacity to within one 
bit as Ps larger than 10, and is within a small fraction of a 
bit for Ps greater than 100. 

2) Fig. 9 also shows the achievable ANC rate by setting 
the amplification gains to the upper bounds of the power con- 
straints given in Lemma 1. There is a constant gap approximate 
1.5 bits from the upper bound. 

Example 2: 4-Layer Network with n relays at layer 2 

We next present the performance of ANC in the 4-layer 
network with n relays on layer 2 depicted in Fig. 10. As the 
number of the relay nodes n sufficiently large, the generalized 
high-SNR regime can be satisfied. 

The result shown in Fig. 11 implies that the gap between 
the achievable ANC rate and the upper bound decreases within 
one bit as n lager than 5, and within a small fraction of a bit 
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Fig. 12. Gaussian two-hop multiple access channel 



network is given as follows. 
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To distinguish a specific AF scheme B, we use {B} to indicate 
Fig. 11. Upper bound of ANC capacity and increase of achievable ANC all feasible relay amplification gains are allowed. 

rate with number of relays at layer 2 in 4-layer network. t ^i • -i ^i • j • ^ ^ ^i j ^..■ 

In the similar way, the received signal at the destmation 

can be represented by an output of an equivalent scaler MAC 

shown as follows, 
for n greater than 40 as the transmitting powers of the nodes 

on layer 2 are upper bounded by 2. 



TJD = XSi,eq + XS2,eq + Zeq, 



(50) 



VI. Extension to Two-hop MAC Network 

In this section, we progressively study the performance of 
the AF scheme in a two-hop parallel MAC depicted in Fig. 12. 
The structure of the optimal rate region of this network is first 
characterized by an asymptotical approach. Then a practical 
low complexity AF-based ANC scheme is proposed. 



A. Achievable Rate Region of Two-Hop MAC via AF 

For simplicity, we focus on a two-source case. With the help 
of the relay nodes, the two source nodes communicate with the 
destination node individually. Assume each non-source node 
introduces an i.i.d. Gaussian noise drawn according to the 
distribution Af(0,l). With the superposition property of the 
wireless channel, the received signals at each relay node is the 
linear combination of the signals from the two source nodes 
and the Gaussian noise. To perform ANC, the relay nodes 
adopt the AF relay scheme. Hence the received signal at the 
destination node is given as follows. 

yo = hf Bhoixs, + hf Bho2a;s2 + hf Bzi + zd, (48) 

where xsi and xs2 are the transmitting signals from the two 
source nodes, hoi and ho2 denote the broadcast channels 
between the two source nodes and relay nodes respectively, 
and hi represents the multiple-access channel between the 
relay nodes and the destination node. The other parameters 
are defined in sections II and III. Again, assume each relay 
node has an individual transmitting power constraint (2). As in 
the point-to-point case, it is easy to verify that the equivalent 
constraint of the power constraint at each relay node k in this 



where xs^^eq = hfBhoix^i and xs^.eq = hfBho2a;52 are 
two equivalent transmitting signals of the two sources respec- 
tively, and Zeq = hJUzi+ZD is the equivalent Gaussian noise. 
With a specific B, the two-hop MAC can be considered as a 
conventional Gaussian MAC. To distinguish them, we denote 
the former one as MAC(B). As the well-known result of Gaus- 
sian MAC, the source nodes adopt the Gaussian codebooks. 
The r.v.'s used to generate the codebooks are xs^ ^ A/'[0, PsJ 
and xs2 ^ J^[0,Ps2]- The codebooks consist of [2"^^] and 
|"2ni?2"| codewords of length n respectively, and the decoding 
error probability tends to zero as n — > oo. 

Corresponding to different AF schemes employed by the 
relay nodes, MAC(B) may have different achievable rate 
sets denoted by 7?.(B) as shown in (51), where A — 
Ilndiag{Ps,,Ps2}'tio^ and Hq = [hoi,ho2]. Denote by 
TZ ({B}) the union of all achievable rate sets. However, note 
that the union region may not be convex. Therefore, to obtain 
the entire achievable rate region, we consider the time-sharing 
technique according to different AF schemes. Consequently, 
the following corollary is established. 

Corollary 1: The optimal achievable rate region of the two- 
hop MAC via AF relay scheme is the closure of the convex 
hull of all the possible rate sets (51), 

col(coIlv(7^({B}))). 

The convex hull is taken by the time-sharing technique 
among all the rate sets. However, it is even more complicated 
to obtain the optimal achievable rate region of the two-hop 
MAC explicitly than to obtain the optimal rate as in the 
point-to-point case. So we focus on giving an asymptotical 
characterization of the rete region in the sequel. 



n{B) ^ \ (Ri, R2) : Ri < C 



hfBhpihg^iBhiPsi 
hI'BBhi + 1 



R2<C 



hl^Bhoah^aBhiPs. 



,Ri+R2<C 



f hI'BABhi 



7^^*(B)-Ui^l,i?2):i?l<C 



hf BBhi + 1 y ' -"i ' ^"-^ - - Vh^^BBhi + 1 
hfBhoihJiBhiPs, 



h:[^BBhi 



i?2<C 



hfBho2^2BhiPs. 



,Ri + R2<C 
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hfBBhi 



hfBBhi 

7^r* (B) = {(i?i,i?2) : i?i < C (hfBhoihJiBhiPsO , 

P2 < C (hf Bho2hj2BhiPs,) , i?i + P2 < C (hf BABhi) } 



(51) 



(52) 



(53) 



B. Outer Bound to Achievable Rate Region 

To investigate the structure of the optimal rate region, we 
first derive an outer bound for it in the similar manner as 
we obtain the upper bound of the ANC rate for a layered 
relay network. If omitting the received noise at the destination 
node or the propagated noises from the relay nodes, two larger 
rate sets T^J"* (B) and T?.™* (B) can also be obtained with 
AF scheme B respectively as shown in (52) and (53), which 
yield an outer bound for the optimal rate region given in the 
following corollary. 

Corollary 2: An outer bound of the achievable rate region 
of the two-hop MAC via AF is the intersection of two outer 
bounds shown as follows. 

col (conv (Ter* ({B}))) fl col (conv (7^^"* ({B}))). 

In the following, we wish to address the issue on de- 
riving the two outer bounds, i.e., col (conv (7?.°"* ({B}))) 
and col (conv (7?.2"* ({B}))), in the closed forms. The first 
observation is that the second outer bound can be easily 
expressed in a closed form. It is easy to verify that the three 
Gaussian capacities in the rate set 7^2"* (B) can be maximized 
simultaneously by a specific choice 



B* 



diagi^lS^'^ ' 



■P'/' 
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Up\ 



(54) 



The corresponding rate set is obtained, which covers any other 
rate sets in W^"^* ({B}). Moreover, note that since 7?.2"*(B*) 
itself is a convex set, the second outer bound can be given 
in (55). However, to obtain the first outer bound is not such 
straightforward. Several useful lemmas are given as follows 
for characterizing it. 

Lemma 2: To have the first outer bound of the achievable 
rate region of a two-hop MAC, it is sufficient to take x = Bhi, 
which is called the amplification vector hereafter, in the linear 
subspace span {hgi, « = 1, 2}. 

Proof: The proof is given in appendix D. ■ 



As shown in (52), the three Gaussian capacities are inde- 
pendent of the norm of x. By absorbing the norm of hoi and 
ho2 into Psi and Pj^ respectively, it is assumed without loss 
of generality, that hoi, i — 1,2 and x are normalized vectors. 
Specifically, a two-hop MAC with two relay nodes is analyzed 
first. By Lemma 2, it is sufficient for us to work on a two- 
dimensional plane. We rewrite hoi, i = 1, 2 and x as 



hoi 



cos a, smo; 



ho2 = [cos/3, sin/3] 
x = [cos 6*, sin 6] 



(56) 
(57) 
(58) 



The key insight of the above expressions is that we can 
find a set of orthogonal bases of space span {hoi, i = 1,2}, 
e.g., Ui = [1, 0] and U2 = [0, 1] , then project all the three 
vectors onto them. Therefore, from the geometric aspect, the 
elements in each vector can be interpreted as the coodinates 
under such base. Assume without loss of generality that a < j3 
in the sequal. By substituting (56)-(58) into (52), the rate set 
7?.J"*(B) can be recast as 



{(Pi,i?2) ■■Ri<c{cf>i{e)), 

R2<C{cj)2{9)),Ri 



i?2<C(0(0))}, (59) 



where 6* e [0,27r], 

01 (0) ^ Ps, cos2 



•-a), 



02 (0) - Ps, cos2 {9-/3), 

(j, {6) = Ps, cos^ {e-a)+ Ps^ cos^ (6* - /3) . 

To emphasize this, later we use 7^°"* (9) to denote the rate 
set in the first outer bound and 7?.°"* {{9}) to denote the 
union of them. With the assistant of the above expression, 
we investigate the relationship between different rate sets in 
the first outer bound and find an interesting result with which 
the scope of 9 to be considered is narrowed to [a,/3]. 



col(conv(7e°"*({B}))) ={(Pi,P2) :-Ri <C(hfB*hoih^iB*hiPsJ , 

P2 < C (h[B*ho2hj2B*hiPs,) , Pi + P2 < C (hf B* AB*hi) } (55) 



Lemma 3: To have the first outer bound of the achievable 
rate region of a two-hop MAC with two relay nodes, it is 
sufficient to take 9 E [a,/?]- 

Proof: The proof is given in appendix E. ■ 

Therefore, we only consider 7?.°"* {{&}) for 9 € [a, /3] in the 
sequel. There still another problem remains to be solved. Is the 
time-sharing operation necessary? In other words, is the union 
region is a convex one? From (59), it shows that if both (j)i (9), 
i — 1,2 are concave, then (f) (9) is concave, and since C (•) is an 
increasing concave function, the three capacities in (59) are all 
concave functions of 9. Consequently, the union of all rate sets 
"^r* ({^}) is a convex set. In this case, conv (7^^* {{0})) = 
jiout ({^1). In the following lemma, a sufficient condition for 
this case is proposed. 

Lemma 4: The first outer bound of the achievable rate 
region of a two-hop MAC with two relay nodes is convex 
when \a- f3\< f. 

Proof: By Lemma 3, we only consider 9 e [a, /?]. Hence 
if \a- (i\< f , then \9~a\<^ and |6i - ^| < f . From (59) 
we see that 



^01 {9) = -2Ps, cos 2 (0 - a) < 



(60) 



(61) 



when |6' — q;| < j, which implies (t)i (9) is concave, and 

^h{0) = -2Ps,co82{9-p)<0 

when |6' — /3| < j, which implies 02 (9) is concave. 

Consequently, (p (^) is concave, and thus 7?.°"* ({^}) is 
convex. ■ 

Indeed, each result obtained in the above lemmas can be 
extended to a two-hop MAC with n relay nodes directly. The 
key step is the observed in proof of the following theorem. 

Theorem 5: To have the first outer bound of the achievable 
rate region of a two-hop MAC with arbitrary n relay nodes, it 
is sufficient to take 9 <E [a,/?]. If | arccos(hj^iho2)| < f , then 
jiout (1 5/}) is a convex set and hence 



conv (7^^* {{9})) = Ter* ({^}) , ^ e [a, (3] 



(62) 



Proof: The proof is given in appendix F ■ 

However, the structure of the first outer bound is still 
unknown in general case. We progressively work towards 
this goal. A commonly used method to characterize different 
rate-tuples on the boundary of a multiuser capacity region 
is via solving a sequence of weighted sum rate maximiza- 
tion problems, each for a different nonnegative rate weight 
vector of users (see [14] and also [21]). For a two-hop 
MAC, such technique can slao be employed to character- 
ize TZi"* {{&})■ That is to solve an optimization problem 
max^i_Ri + /i2-R2,MiiM2 > 0. It is clear that the rate pair 
{Ri, R2) on the boundary of the union region must be on the 
boundary of some rate region 7?. J"* (6*), which should satisfy 
(52). By such observation, we first analyze the relationship 
between the boundary points of 7?.°"* (6') and the given weight 
vector {111,1x2) in the optimization problem. Without loss of 
generality, assume that /ii and 1x2 are normalized in [0, 1], and 
that /ii + ^2 = 1- Therefore, the maximization can be recast 
as max^i?i +/2i?2, fi = fii,jl = ^2- For any given 9 e [a, P], 
it is easy to verify that 



upper-diagonal corner point 




lower-diagonal comer point 



Fig. 13. Rate set in first outer bound 

• case 1: when /i ~ 1/2, the rate pairs maximize the 
sum rate in rate set 7^°"* (^) ^re the solutions of the 
optimization problem, which satisfy 

Ri{9) + R2{9)^C{(l){9)). 

• case 2: when ^ = 1, the rate pairs maximize the 
individual rate i?i in rate set 7?.°"* (9) are the solutions 
of the optimization problem, which satisfy 

Ri{9) = C {cf>i (9)), < R2{9) <C{<P (9)) - C (0i {9)). 

• case 3: when ^ = 0, the rate pairs maximize the 
individual rate R2 in rate set 7?,°"* (9) are the solutions 
of the optimization problem, which satisfy 

< Ri{9) <C{<P {9)) - C (02 [0)), R2{0) - C (02 {e)). 

• case 4: when ^ e (O, ^), the rate pair at the "upper- 
diagonal" corner point of the pentagon region (see Fig. 
13) in rate set 7?.°"* (9) is the solution of the optimization 
problem, i.e., 

Ri{9) ^C{<P [9)) - C (02 [0)), R2{e) - C (02 [9)). 

• case 5: when /i G (5,1), the rate pair at the "lower- 
diagonal" corner point of the pentagon region (see Fig. 
13) in rate set 7?.°"* {9) is the solution of the optimization 
problem, i.e., 

Ri{9) =. C (01 {9)), R2{9) ^C{cj) (9)) - C (0i (9)). 

With the analysis made above, we conclude that for any 
given /x, the boundary points of 7^°"* {{9}) can be derived 
by miximizing nRi{9) + flR2{9) with respect to 9 and the 
optimal solution is denoted by 9*{fi). We explicitly obtain 
the rate pairs on the boundary which maximize the sum rate 
and the individual rates of 7?.°"* {{9}) in the sequel. Then 
two equations will be presented, whose solutions maximize 
liRi{9) + p,R2{9) for ^ E (O, i) and /i G (|, l) respectively. 

Let us start with the maximum sum rate of 7?. J"* ({^})- The 
following lemma is established. 

Lemma 5: The maximum value of i?i + R2 of 7?.J"* ({^}) 
is given by C (0 {()*{^))), where 9*{^) is shown in (63) on 
the top of the next page. 

Proof: For any 9, the maximum sum rate R1+R2 is upper 
bounded by C(0(6')). It corresponds to the maximization 
problem when /^ = 1/2 as shown in case 1. By setting the 
derivative 0' (9) to zero, we have 

0' (6*) = -Ps, sin 2 (9 -a)- Ps^ sin 2 (0 - /3) = 0. (64) 

Then, it follows that the solution of the above equation satisfies 

+ 0/3*^1 \ ^Si sin2Q: + Ps2sin2/3 

tan 26^ (-) = • (65) 

^2' Psi cos 2a + Ps2 cos 2/3 



,,1, 1 /Ps, sin 2a + Ps, sin 2^ 

^2^ 2 \Ps,cos2a + Ps; cos2p 



-'h^-\ 



arctan 



Ps^ sin 2a + Ps^ sin 2/3 
Psi cos 2a + Ps2 cos 2/3 



Psi sin 2q: + Ps^ sin 2^ 
Pgj cos 2q: + Ps2 cos 2/3 
Psi sin 2a + Ps^ sin 2/3 
Pgj cos 2a + P52 cos 2/3 



>0 



<0 



(63) 



ti [Ps^Ps2 sin 2 (/3 - e) cos^ [6 - a) + Pg, sin 2 (0 - a) (l + Pg, cos^ [0 - p))] 
= fi [Ps,Ps2 sin 2 (/3 ~ 61) cos^ {0 - a) + Ps^ sin 2 (/3 - 6I) (l + Pg, cos^ {0 - p))] 

M [PsiPs2 sin 2 (61 - a) cos^ (6' - ^) + Ps^ sin 2 (6* - a) (l + P^^ cos^ {0 ~ a))] 
= M [PsiPs2 sin 2 (61 - a) cos^ {9 - l3) + Ps^ sin2 (;3 - ^) (l + Psi cos^ {9 - a))] 



(68) 



(69) 



Since 0*{^) G [a,/3], 6'*(i) is obtained as shown in (63) on 
the top of the page. The AF scheme corresponging to 0*{^) 
is denoted by B^^'^'. By the notations used in the proof of 
Theorem 5, 



B(io) 



.(10) 



diag 



H7 



,1 



1, 



cos6'*(-)ui +sin0*(-)u2 



where the constant c'^"' is chosen such that the constraint (49) 
is satisfied and at least one of the amplification gain achieves 
the upper bound. ■ 

Then we consider the maximum individual rates of 
■j^out ({5»}). The following lemma is established. 

Lemma 6: The maximum value of individual rates Ri, i = 
1, 2 of 7^^* {{0}) are given by C {(f)i (a)) and C {(j>2 (/3)). 

Proof: To derive the maximum individual rates is equiva- 
lent to solving the maximization problems maxC((/!)i {0)) and 
maxC((/!)2 {0)) according to /i = 1 and ji — Q respectively as 
shown in case 2 and 3. By setting the derivative 4)\ [0) to zero, 
we have 

(/)'i {0) = -Ps, sin 2{0-a)=O. (66) 



The solution of the above equation is 0*(1) = a. 
By setting the derivative (f)'2 (0) to zero, we have 



<f'2{()) 



Ps2 sin 2 I 



/3) = 0. 



(67) 



The solution of the above equation is 0*{Q) = /3. 

By the notations used in the proof of Theorem 5, the AF 
schemes corresponging to 0*{1) and 0*{O) are denoted by 
B(ii) ^ c^^^) diag [B.^\cos0*{l)ui+sm0*{l)u2)], and 
B(i2) ^ c(i2)(iiag [Hj^i (cos6'*(0)ui + sinr (0)u2)], where 
the constant c'^^'^\i = 1,2 is chosen such that the constraint 
(49) is satisfied and at least one of the amplification gain in 
B^^'^ achieves the upper bound. ■ 

Lemma 7: The weighted sum rate iiRi{0) + J2R2{0) is 
maximized by 0* (/z) satisfying (68) given on the top of the 
page for /i S (O, ^), and is maximized by 0*{ii) satisfying 
(69) given on the top of the page, for /i e (^j l)- 

Proof: Since for /z S (O75)' the rate pair (^1,^2) on 
the boundary can be determined by solving the maximization 
problem /i (C ((/) (6*)) - C {(f>2 (9))) + flC {(f)2 (9)). By setting 
the derivative of the object function to zero, we obtain (68). 




</;<! 



Fig. 14. First outer bound 

By symmetry, the proof of the second part is exactly the 
same as the first part. ■ 

From (68), it is clear that each 9*{^), ^ e (0,^) is in 
[6'*(i),/3] and conversely for each 6* e [0*(i),/3] there exists 
a /i € (0, i) satisfying (68). Consequently, the boundary points 
determined by the optimization problem with ji E (0, 5) 
(see Fig. 14, be) can also be determined by the parametric 
function of the "upper-diagonal" corner point of the region 
corresponding to the rate set Tll^\9) with e [6'*(i),/3]. 
Similarly as the previous case, it is clear that each 0* (/i) for 
given /i e (5,1) is in [a, 0*(i)] and conversely for each 
e [a,6'*(|)] there exists a ^ e (i,l) satisfying (69). 
Therefore, we claim that the boundary points determined by 
the optimization problem with ji £ (5,1) (see Fig. 14, cd) 
can be determined by the parametric function of the "lower- 
diagonal" comer point of the region corresponding to the rate 
set7^^*(^) with e [a,6'*(i)]. 

Finally, the following lemma is proposed for the purpose of 
analyzing the convexity of 7^°"* ({^D- 

Lemma 8: The parametric function of the "upper-diagonal" 
corner point is a concave function of 0, for E [6** (i), /3] , and 
the parametric function of the "lower-diagonal" comer point 
is a concave function of 0, for E [a, ^* (^)] . 

Proof: The proof is presented in Appendix G. ■ 

By the lemmas at hand, the entire region of the first outer 
bound can be fully characterized and thus the outer bound. 

C. Dynamic AF-Based ANC Scheme 

We focus on designing a practical AF-based ANC scheme 
in the left of the paper. Consider the three Gaussian capacity 
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formulas in the achievable rate set 7?.(B), we find that each 
of them corresponds to a point-to-point two-hop relay net- 
work discussed in section V. According to different scenarios, 
differently many AF schemes are stored in the relay nodes 
constituting a dynamic ANC scheme. By employing the time- 
sharing technique, the achievable rate region of the dynamic 
ANC scheme is the closure of the convex hull of the rate 
sets it contains. It is obvious that one can obtain a larger 
rate region by storing more AF schemes in the relay nodes, 
but it will increase the complexity of the relaying scheme. 
So the practical ANC design involves a trade-off between the 
performance and complexity. 

In the following, we consider four typical AF schemes. 
The first three ones are B^'^"\ B(") and B^^^). The forth 
AF scheme considered, denoted by B*^^', consists of the 
amplification gains as the upper bounds of the constraints (49). 
Similarly, the upper layer noises under different AF schemes 
are denoted by Wj^^ , i ^ 0,1,2 and Wj^'. We redefine the 
domination of the noises received at the destination. 

Definition 6: We call the upper layer noise dominates the 
total noises received at the destination in a two-hop relay 
network with AF relays when 



E 



i^'^f 



> E 



(l«)^2 



(Wd ) 



> l,i = 0,l,2. 



Reversely, we call the lower layer noise dominates the total 
noises when 
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Fig. 15. Two-hop MAC with two relay nodes 
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(70) 



> 1 



where (a) follows from the assumption that E 
and (b) follows from Lemma 5 that the maximum sum rate of 
the optimal achievable rate region is no larger than the one of 
the first outer bound, i.e.. 



i>E {w^s^r 



>E {w^j^'^f],i^0,l,2. 



{R1+R2 



Under these two special scenarios, the performence of two 
dynamic ANC schemes are analyzed. 

Theorem 6: For a two-user two-hop MAC with AF re- 
lays under individual power constraints, the gap between the 
achievable rate region and the optimal AF rate region is at 
most 1/2 bit with respect to the maximum individual rates 
and the sum rate when either the upper layer noise or the 
lower layer noise dominates the total noises. 

Proof: For the first scenario, we compute the achieveble 
rate region with a dynamic ANC scheme consisting of three 
schemes, i.e., B^^^^ B*^^^) and B^^^^. The corresponding 
achievable rate region is the closure of the convex hull of 
7?,(B*^^*^ ), i = 0, 1, 2. For the gap between the outer bound and 
the maximum achievable individual rates, the proof is exactly 
the same as the point-to-point case given in Theorem 3. For 
the sake of brevity, we omit the detailed proof here. 

Denote the maximum sum rate of the optimal AF region by 
[Ri + R2)max- We need to show that the maximum sum rate 
of the achievable rate region is at most 1/2 bit from (i?i + 
R2)max- From (51), the sum rate corresponding to scheme 
B(i2) is given as follows. 



-Rl + i?2 
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h:fB(i2)B(i2)hi + 1 

hfB(i2)AB(i2)hi/^[( 
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<c(^(r(-))) 
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■log 1 
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For the second scenario, the dynamic ANC scheme only 
consists of the forth scheme B^^^ The proof is almost the 
same as the previous case. Therefore, the proof is skipped 
here. ■ 

To be more general, we can obtain a more complicated ANC 
scheme which is not limited to consisting of the above four 
AF schemes. It is clear, the more AF schemes stored at the 
relay nodes the lager the achievable rate region we may get. 
Hence, according to different application goals, we can have 
a tradeoff between the complexity of the coding scheme and 
the performance of it as the common argument. 

D. Examples 

A two-hop MAC with two relay nodes is depicted in Fig. 14. 
Consider the general case and the two cases either the upper 
layer noise or the lower layer noise dominates the total noises 
received at the destination node. The corresponding network 
parameters are given in Figs. 15-17. 

In the first case, we employ an ANC scheme with the four 
AF schemes B'^*',i — 0,1,2 and B^^' and the time-sharing 
technique. We explicitly plot the four rate regions in Fig. 15. 
The numerical result shows that the gap between the inner and 
outer bounds is approximate to 1/2 bit. 

In the second case, we employ an ANC scheme with the 
four AF schemes B'^^*\« = 0,1,2 and the time-sharing 
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Comparison of two outer bounds and time-sharing inner bound 



Comparison of two outer bounds and inner bound 
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Fig, 16, Time-sharing inner bound with AF schemes 1-4 and outer bound Fig, 18. Inner bound with AF scheme 4 and outer bound 
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Fig, 17. Time-sharing inner bound with AF schemes 1-3 and outer bound 



technique. Since E 



w 



(2)^,2 
D I 



E 



(ll)^2 



«'') 



= E 



/ (12)\2' 



E 



i< 



(10)n2 



6.33, and 



= 3.29, the upper layer noises 
doihinate the total noises. Therefore the inner bound obtained 
by the ANC scheme yields a 1/2 bit gap to the outer bound 
for the sum rate and two individual rates as shown in Fig. 16. 
In the last case, only the forth AF scheme is em- 



ployed. Since E 



E 



{w 



(11) 

D 



(2)n 



E 



(w 

{^D ) 



E 



{< 



(10)n2 



0.38, and 



0.20, the lower layer noise 
dominates the total noises. Therefore, the inner bound obtained 
by the ANC scheme yields a 1/2 bit gap to the outer bound 
for the sum rate and two individual rates as shown in Fig. 17. 

VII. Conclusion 

We cast the problem of computing the maximum achievable 
rate of AF-based ANC scheme as an optimization prob- 
lem. For a multi-hop wireless relay network, assuming an 



individual power constraint at each relay node, we derive 
a tighter amplify-and-forward upper bound for the optimal 
transmission rate. The closed-form of such bound is derived 
for a layered network with a single source-destination pair. 
We have provided the analysis of the AF-based ANC scheme 
both for a point-to-point case and a two-hop MAC case with 
the assumption of instantaneous relay operation. However, in 
general wireless relay networks with AF relaying, the resulting 
input-output channel between the source and the destination is 
an intersymbol interference (ISI) channel with colored noise. 
We want to extend the results in this paper to that general 
case, and we conjecture the corresponding ANC scheme may 
be much more similar to the conventional network coding 
framework. 

Appendix A 
Proof of 'L^mmk 1 

Proof: Let Wk be the total noise received at node k. 

where {ifc} denotes the nodes appearing in the relaying path 
from S to k. The signal received at node fc is expressed as 

Vk = hM^s + Wk (72) 

We prove the theorem by induction. Consider first the node k 
whose neighboring nodes set only contains the source node S. 
From (1) and (3), 



(a) 



E [xl] = f^lE [y2] < 



^k 



{l + 5k)PR,k 



{h%kPs 



W pc/p 



l)=P, 



(73) 

where (a) is obtained from the condition (11), and (b) is drawn 
from (9). To prove that the theorem holds for any node, we 
assume that all the neighboring nodes of node k satisfy the 
power constraints such that 



E [x]] - E /3^ ifs^.xs + W,y 



<Pf^ieV(fc) (74) 
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Then we consider the transmitting power at node k with the 
amphfication gain chosen satisfying the condition (11), 



It concludes that the maximum of SNRi^^ is no less than the 
maximum of SNR. We first point out that 



E[4]'^E 



(b) 



Pi I fs.kxs + ^ Pjhj^kWj + Zk 
jev(k) 



7-B,„g-GA„, 



(77) 



where 



< coo < fs.kPs + E 



uev(fc) 



l„) = [gi---gj---9n,A : 



^'^coo<(coi+ E hlkP^,[fl,Ps + E{w])]+l 



jevik) 



id) 



<coo<coi+ J2 ^-fe-Pf^ + l 



jevik) 



(e) 



<coo< E Piho,kPihi.k ^E{w^)E{wf) 



j,jev{k),i^:j 

~fs,fsAPs]+ E ^-fc^'^' + l 
jev(fc) 



Bio = diag {/3i ■■■fiy /3„,^ } , j G £;, 
g=(hLiBL_i---H 

G = diag{gi---5j •••g„,J. 
The received signal vector at layer Zq is 

y(o = [/s,i • • • /sj- • • • /s,n,J a;s + Z;^ 
and denote by 

P(o = [/s,i • • • fs,j ■ ■ ■ /s,n,J yPs 

A solution to 

max SNRig 



(78) 



(79) 



(80) 



(/) 



can be found by representing the object function in the form 
of Rayleigh quotient that would then be maximized. With the 
notations defined above, SNRi^^ can be expressed as 



< ^00 E ''hp^' + E hkKk JpT'pr' + 1 



jev(fe) 



1^3 



Up 



p; 



where cqo = 



(75) 



P, 



Up 



SNRi, (7) = 
The optimal solution is 

lopt 






7^7 



(81) 



(82) 



Coi = E I^J^J,kl3lhLk [E {WjWi) + fs.jfs.jPs], 



where c is a constant which is chosen such that the power 
constraints are satisfied, the corresponding optimal value is 



!,jev(fc) 



5iVi?io(7op*)=Pi^Pio, 



(83) 



and 



(a) follows from (6), (7)and (11), 

(b) follows from (10), 

(c) follows from (6), 

(d) follows from the assumption (13), 

(e) follows from the Schwarz inequality, 

(f) also follows from (13), and the fact that 

2 



E{w])Eiwf) + fs,,fsjPs 

< [E H) + fhPs] [E H) + fhPs] 
Then we complete the proof. 



Appendix B 
Proof of Theorem 1 

Proof: We denote the SNR function at the destination 
of an ideal network where only the nodes at layer ^o^ 'o = 
1,2 ■ ■ ■ L, introduce the Gaussian noises as SNRi^^. With the 
identical ANC scheme, it is obvious that 



which can be interpreted as the total received signal power at 
nodes in layer Iq. 

It follows that the larger the elements in Pjj,, the larger 
the optimal value. From (9) in definition 3, we know that the 
received signal power at node k in layer Iq is upper bounded 
by Prm- Therefore the total received signal power at layer 
lo is upper bounded by P|J ; P_r,;o- Consequently, the upper 
bound of SNRi„ (lopt) can be given. 



SNRrfh. 



opt) 



^RAo^R^loi 



(84) 



and the corresponding Gaussian channel capacity evaluates to 



R 



Up 



1 



^Og{l + PliPR,lo), 



(85) 



which can be served as an upper bound to the optimal achiev- 
able rate of the original network. By taking Iq — 1,2- ■■ L, a 
set of upper bounds can be derived which leads to the result 



1 



R^P = min rY,^ = _ min ^ ^ log (l + PIi„Prj„) 



lo = l,2---L 



la = l.2---L 2 



SNRi„ > SNR. 



such that Ropt < R^^ holds. 



(86) 



(76) Then we complete the proof. 
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Appendix C 
Proof of Theorem 2 

Proof: It is obvious that the maximum achievable rate 
Ropt is no less than the rate derived from the scheme proposed 
in this paper. To present the results, C{x) is used to denote the 
capacity of Gaussian channel with SNR x. 



( 



Ropt ^ C 



E 



(h];_iBL-i---HiBihoa;, 



\ 



E 



(a) 
= C 



1=1 

\ 



E[zl] 



cf [Pkk^RM, 



(1 + So) 



lo~l 



C4 + cfPij,„Pfl,io + l 



V 



/ 



(b) 

>c 



^R,lo^R,lo 





(l + <^o)'""' 


1 - 


1 

{i+5oy'>-\ 


Pl^P R,1,+CS 



(87) 



where C4 — E 



. 1=1 
We first claim that (a) follows from 



+ C2. 



E [(hLiBL-i • • • H,„B,„z,„)'J = E [g^B,„z,„z;^B,„g] 

= E [g^B,„z,„z^B,„g] = A^GGAo = cIpI^^Prj,, 

where g and G are defined in (77) and the last equality follows 
from (34), and 

^[(hLiBL-i---HiBihox,)2] 

(1 + 6ot'-' 
Then, we claim that (6) follows from 



R,loi 



E 



/io-l 



J^i^l-i'^L-i-'-HiBiZi 
1=1 

Y,E[{hl_,BL-i---liiB 



\ 1=1 

lo-l 



m 



1=1 

In-l 



< 



^—1 \ J- 



<5o 



lo-l 



^i^L-l^L-l ■ ■ -HioBioPr.I 



E 



-50 



-ACPlio^RJo) 



=1 (1 + -^o)' 

where the inequality follows from Lemma 1 in [15]. 

By taking Zq = Ij 2 • • • L, a set of lower bounds can be 
derived which leads to the result. ■ 



Appendix D 
Proof of Lemma 2 

Proof: To prove Lemma 2, we justify the following 
proposition. To have each rate set 7^°"*(B), it is sufficient 
to take X = Bhi in the linear subspace span {hoi, « = 1, 2}. 
Then Lemma 2 can be viewed as an directed corollary of such 
proposition. 

On one hand, any n-dimensional vector x can be decom- 
posed as 

x = xi+X2, (88) 

where xi e span {hoi,i ~ 1, 2}, X2 G span^ {hoi, « = 1, 2}, 
and thus x^X2 = 0. 

On the other hand, span-'- {hoi,i — 1,2} is the solution 
space of the homogenous linear equations Hj^x = 0. There- 
fore, it follows that hj^jX2 = 0, i = 1, 2. 

Therefore, from (52), it is easy to verify that given any 
amplification vector x = Bhi, there exists another one y = 
B'hi = cxi such that 7e5'"*(B) C 7^^*(B')■ Note that the 
constant c ensures the scheme B' in the feasible set {B} and 
has no effect on the rate set 7?.°"*(B'). 

In other words, the corresponding vector component of the 
amplification vector belonging to span^ {hoi, J = 1,2} will 
not contribute to the received SNR for the three Gaussian 
capacities in the rate set (52). 

Then we complete the proof. ■ 

Appendix E 
Proof of Lemma 3 

Proof: The lemma can be proved by a consequence of 
careful calculation. Several basic conditions are described first. 
As assumed in the previous section that all the channel gains 
are positive values, the angles a and /3 are both in [O,^]. 
Without loss of generality, assume < a < j3 < ^ hence < 
j3—a < J. From (59), we can see that the outer bound 7?.™* (6) 
only dependents on the inner product of the channel vectors 
and the amplification vector. Although 9 can take all values 
in [— 7r,7r] , it is easy to find that x-^hoi = cos^ {9 ~ a) — 
cos^ {6 — aztn) and the same argument applies to the inner 
product between x and ho2- This implies that the value of 9 
can be limited to [— f , f ] • Then the proof can be completed 
by the following discussions. 
. Case 1. When 9 e [/3, f]. 



X hoi 



cos 



(a) 

-a) < cos^ (/3 - a) , (89) 



where (a) follows hy 0</3 — a<9 — a<^, and 

x^ho2 == cos^ (9- (3) < cos^ (^ - P) = 1. (90) 

Therefore, by taking 9* = (3, n'^^'>{9) C 7er*(^*) holds 
for all 9 e [/3, f ] . 

Case 2. When 9 E [a, (3], there is nothing to prove. 
Case 3. When 9 e [/3- §,a]. 



x^hoi = cos^ (9 — a) < cos^ (a — a) = 1, 



and 



x^hn 



(a) 

cos^ {(3-9) < cos^ {(3 -a). 



(91) 



(92) 
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where (a) follows byO</3-a</3-6'<f. 

Therefore, by taking 9* ^ a, 7^5'"*(6') C 7^5'"*(6'*) holds 

for all 0e [/3- f,a]. 

Case 4. When 6* e [a- f,/3- f], 

let A = /3 - I - 61 e [0, /3 - a] and 61' = /3 - 



Ae 



X' ho2 = COS {P -9) = cos ( — 



2 C^ 

= cos (-- 
1). If 9' < a, then 
x"^hoi — cos^ (a — 61) 

(«) 2/ 

< cos I a — /3 



A 



cos^ {P - 



(93) 



= cos 1 a - 


-^+|.A) 


-f-A) 


= cos^ (a — 9') , 




(94) 


<a-e' = 


a-/3+f -A< 


<f- 





Q;-,3+f + A = Q;- 

Therefore, in such case, for any 9 there always exists a 

61' e [P- f ,a] that 7^^*(^) C 7^^*(6'') holds. Hence 

from the result of case 3, by taking 9* = a, 7^™* (9) C 

7^5'"*(^*) holds. 

2). If a < 9' < 13, then 



A 



X hoi — cos (a — 0) = cos I a — /3 + 



(a) 



f-- 



< cos^ ^l3-a- — + l^] = cos" (^y - a) 

(95) 

where (a) follows from 0<a-9'^a-P+^-A< 

a~ l3+^ + A^a-9 <^. 

Therefore, in such case, for any 9 there always exists a 

9' e [a, 13] that 7^^*(^) C 7^^*(^') holds. 

Cases. When 6* e [-f,a-f], 



x^hn 



= cos^ (6* - a) : 

(a) 2 

< cos (/3 — a) . 



cos 



(96) 



where (a) follows from f-/3< f-a< 9-a + Tr< f, 
and 



X^ho2 



cos^ (6' - /3) < cos^ (/3 - /3) = 1. 



Therefore, by taking 9* = /3, 11°^"* (9) C 7^!j'^ 
for all 0e [-f,a- f]- 
Then we complete the proof. 



(97) 
') holds 



(P^. cos2(0-a),P,^ sm2{0-a)'j 
C 




P,, cos2{ff-a) + P^_ cos2(6i-/?y 
P, sin2(ff-a) + P,, sin2(6'-/?) 



2(e-y9),/>s^sin2(e-y9)) 



Fig. 19. Complex plane representation 1 



(Pj^ cos 2/?,-/',, sin 2/?) 




IP^. cos2Qr,-P^. sin2«) 



( P,. cos 2a + P,, cos 2/9, -P, sin 2a - P,, sin 2/?) 
Fig. 20. Complex plane representation 2 

hoi = cosaui +sinaU2, 
ho2 — cos/?Ui + sin/?U2, 



(99) 
(100) 



where cos 6* = x-^Ui, cos a = hj^ui and cos/3 — h|5"2Ui, and 
a,/3 G [O, ^]. Since all the channel gains are positive, it is 
clear that |a - /3| < |. 

Again, the rate set 7?.(B) for a two-hop MAC with n relay 
nodes can be recast as (59) with (98)-(100). Therefore, by 
Lemma 3 it follows that to have 71°"* {{9}), it is sufficient 
to take 9 E [a, /3], and by Lemma 4 it follows that when 
la- I3\< f , 7^^* {{0}) is convex. 

Then we complete the proof. ■ 



Appendix G 
Proof of Lemma 8 

Proof: The function of the "upper-diagonal" corner point 
of each rate set 7^°"*(^) is given as follows. We consider the 
boundary points determined by 9 G [9*{^), (3]. 



a;-iln(l + 0(0))-iln(l + 02W) 



(101) 



Appendix F 
Proof of Theorem 5 

Proof: As proved in Lemma 2, it is assumed with- 
out loss of generality that the amplification vector x e 
spanjhoi, ho2}. Note that in the proof of Lemma 2, there 
is no restriction on the dimension of x. To apply the same 
argument to obtain (56)-(58), it is assumed that x, hoi and 
ho2 are unit vectors. Then we find a set of unit orthogonal 
bases of linear space spanjhoi, ho2}, denoted by (ui,U2) 
such that 

X = cos6'ui + sin6'u2, (98) 



d^y _ 93 {0} (1 + jO)) (1 + 4>2 {0}) - 92 (0) Oi {0} 



where 



92(0) 

[(l + </-(0))0^(0)]^g?(0)' 
4>' {9) (1 + 02 {9)) 



(102) 



^^^^^" (1 + 0(0)) 0^(0) ^' 

92 {0) = 0' (6^) (1 + 02 {9)) - 0'2 {9) (1 + {9)) , 

g:,{9)^^"{9)M0)-^'{0)^2{O). 



92 (6) =^' (6) (1 + 02 (6)) - ^'2 {0) (1 + -^ m 

= [-Ps, sin 2 (6 -a)- Pg, sin 2 {6 - /3)] [l + Pg, cos^ {9 - /3)] 

+ Ps, sin 2 (6* - /3) [1 + Ps, cos^ {6 - a) + Ps, cos^ (6* - p)] 
= - Psi Ps2 cos^ {9 - /3) sin 2 {6 -a)- Pj^ cos^ (6* - /3) sin 2 (6* - /3) - Pg^ sin 2 {6 -a)- Pg^ sin 2 

+ Ps^ sin 2 (6* - /3) + PsiPs2 cos^ {9 - a) sin 2(0-^) + PJ^ cos^ (6* ~ 13) sin 2 (0 - ;3) 



(0-/3) 



- Ps^ Ps2 cos^ (61 - ;3) sin 2 (e - a) - Ps, sin 2 (9 -a) 



Ps^Ps2 cos^ (S - P)sin2\ 



Ps, sm2\ 



Ps, Ps2 cos^ {9 ~ a) sin 2{9 - p) 
Ps, Ps2 cos^ (6* - a) sin 2 (;3 - 61) 



(103) 



53(^)=0"(e)0^(^)-0'(^)0i'(e) 

= [-2Ps, cos 2 (6* - a) - 2Ps, cos 2 (6* - /3)] [-Pg, sin 2{9~ /3)] 

- [-2Ps, cos 2(9- /3)] [-Psi sin 2 (6* - a) - Pg, sin 2 (6* - /3)] 
=2Ps,Ps2 sin 2 (6* - /3) cos 2 (6* - a) + 2P|^ cos 2 (0 - /3) sin 2 {9 - (3) 

- 2Ps, Ps2 cos 2 (6* - /3) sin 2 (6* - a) - 2P4 cos 2 (6* - /3) sin 2 (6* - /3) 
=2Psi Ps2 [sin 2 (0 - /3) cos 2 (6* - a) - cos 2 (6* - /3) sin 2 (6* - a)] 

= -2Ps,Ps,sin2(/3-a). 



(104) 



54 (e) -</-' (6) 0'2 (0) 

= [-Psi sin 2 (6* - a) - Pg^ sin 2 (0 - 13)] [-P^^ sin 2 (6I - /?)] 
= - Ps2 sin 2 (/3 - 61) [P5, sin 2 (61 - a) + P^^ sin 2{9 - (3)] 



= -Ps2sm2{P-9) 



PI 



Pi + 2Ps,Ps2 cos 2i(3~a) sin 2 (6^ - ^o) 



= -a/p: 



P^ + 2Ps,Ps2 cos 2 (/3 - a)Ps^ sin 2 {(3 -9) sin 2 (6* - 0o) 



(105) 



yl = (P51 cos 2 (6* - a) , P5, sin 2 (0 - a)) , 
B - (P5, cos2 (0 - 13) , Ps2 sin2 (0 - /?)) , 
C == (Psi cos 2 (61 -a) + Pg^ cos 2 (61 - /3) , Pg^ sin2 (6I - a) + Pg, sin2 (6I - /3)). 



(107) 



and 



34 (0) = </-' {0) ct^'2 {0) ■ 



To verify the concavity of the parametric function (101), we 
carefully check each term in (102). The computation of 32 (6^) 
is given in (103) on the top of the next page. We claim that 
52 {0) < for 6* e [r(i),/3]. Because 

2(6'-a)e [0,7r] 
^ - Psi Ps2 cos^ (0 - P) sin 2 (0 - a) < 0, 
- Ps^ sin 2 (6* - a) < 0. 

2(/?-0)e[O,^] 

=> - Ps^Ps2 cos^ (61 - a) sin 2 (^ - 6I) < 0. 
which implies that each term in (103) is non-positive. 



The computation of g^ (0) is given on the top of the page 
(104), we claim that 53 (9) is a non-positive constant because 

2{l3-a)e [0, n] ^ 33 (^) = -2Ps,Ps2 sin2 (/3 - a) < 0. 

The computation of 54 (0) is given on the top of the page 
(105), where the equality (a) follows from the discussion 
below. 



Psi sin 2 I 



a) + Ps, sin 2 (0-/3) 



= Im 



p^^^J2ie-a) ^ p^ 



82^ 



Me-H) 



(106) 



We use A, B and C to denote PsieJ2(e-a)^ Pg^ 6^2(9-/3) ^^^j 

-PsiC-'^^^""^ + PssG-'^'^^'^^ respectively as shown in Fig. 19 
and (107). Hence, 
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/ Ps, sin 2a + Ps, sin 2/3 \ 

BjrctrSjii I -— — _ ^ I , 

V Psi COS 2a + Ps^ cos 2/3 / 

/Ps, sin 2a + Ps, sin 2/3' 

+ arctan — - 

\ Psi cos 2a + Ps2 cos 2/3 ^ 



Pg, sin 2a + Ps, sin 2/3 ^^ 
Psi cos 2a + P52 cos 2/3 ~ 
Pg, sin 2a + Ps, sin 2/3 ^^ 
Psi cos 2a + Ps2 cos 2/3 ~ 



(111) 



PsiSin2( 



'- a 



P 



's, sin 2(6*- /3) = I C| sin 99(C), 

(108) 



where |Cr = Pj^ + P^ + 2P5,P5, cos 2 (a - ^) and ip (C) 
denotes the phase of C. Furthermore, to obtain the phase 
(f (C), we first derive the phase of C" — Ce~^^^. In con- 
sequence, ip (C) — if (C) + 26. As shown in Fig. 20, the 
coordinate of C" is given as 

C" = (Psi cos 2a + Ps^ cos 2/3, -Pg^ sin 2a - Ps^ sin 2/3) . 

(109) 
Therefore, 



tan (93 (C)) 



Ps^ sin 2a + Pg^ sin 2/3 



(110) 



Psi cos 2a + Ps, cos 2/3 ' 

hence ip{C') = -2e*{\) (compare (111) and (63)). There- 
fore, we conclude that equality (a) holds. Applying the same 
argument as we justified g2{0) < , for e [0*(i),/3], 

2 [e - r (i)^ e [0, tt] ^ sin 2 (^0 - r (i)^ > 0, 

2(/3-6') e [0,7r] ^sin2(/3-6') > 0, 

therefore, 54(6') < 0. 

Finally, substituting the above results into (102), we obtain 



dx^ 



<0, 



which implies that the parametric function (101) is concave. 

The same argument also applies to the boundary function 
determined by the "lower-diagonal" points for 9 G [a, ^*(|)]. 
Therefore the proof is omitted. 

We complete the proof. ■ 
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